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Abstract. 

Supernovae of type la (SNe la) are very important for cosmography. To exclude 
systematic effects in linking the observed light of distant SNe la to the parameters 
of cosmological models, one has to understand the nature of supernova outbursts 
and to build accurate algorithms for predicting their emission. 

We review the recent progress of modeling the propagation of nuclear flame subject 
to numerous hydrodynamic instabilities inherent to the flame front. The Rayleigh- 
Taylor (RT) instability is the main process governing the corrugation of the front 
on the largest scales, while on the smallest scales the front propagation is controlled 
by the Landau-Darrieus instability. 

Based on several hydrodynamic explosion models, we predict the broad-band UBVI 
and bolometric light curves of SNe la, using our ID-hydro code which models 
multi-group time-dependent non-equilibrium radiative transfer inside SN ejecta. We 
employ our new corrected treatment for line opacity in the expanding medium, 
which is important especially in UV and IR bands. The results are compared with 
the observed light curves. Especially interesting is a recent 3D-deflagration model 
computed at MPA, Garching, by M. Reinecke et al. 



1. Introduction 

Supernovae of type la (SNe la) are important for cosmology due to 
their brightness. They are not standard candles, but can be used for 
measuring distances (i.e., for cosmography) with the help of the peak 
luminosity - decline rate correlation, established by Yu.P. Pskovskii [28] 
and M.M. Phillips [26] (see the review [22]). The physical understanding 
of the Pskovskii-Phillips relation is crucial for estimating the validity of 
cosmological results obtained with SNe. To exclude systematic effects in 
linking the observed light of distant SNe la to the parameters of cosmo- 
logical models, one has to understand the nature of supernova outbursts 
and to build accurate algorithms for predicting their emission. 
This involves: 

— understanding the progenitors of SNe la; 

— the birth of thermonuclear flame and it accelerated propagation 
leading to explosion; 
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— light curve and spectra modeling. 

Although there is no doubt that an SN la outburst is a result of 
thermonuclear explosion, the details of the process are not yet clear. 
We point out some problems which seem most important to us. 

We review the recent progress of modeling the propagation of nuclear 
flame subject to numerous hydrodynamic instabilities inherent to the 
flame front. The Rayleigh-Taylor (RT) instability is the main process 
governing the corrugation of the front on the largest scales, while on 
the smallest scales the front propagation is controlled by the Landau- 
Darrieus instability. 

Based on several hydrodynamic explosion models, we predict the 
broad-band UBVI and bolometric light curves of SNe la, using our ID- 
hydro code which models multi-group time-dependent non-equilibrium 
radiative transfer inside SN ejecta. We employ our new corrected treat- 
ment for line opacity in the expanding medium, which is important 
especially in UV and IR bands. The results are compared with the ob- 
served light curves. It seems that classical ID thermonuclear supernova 
models produce the light curves that fit the observations not so good 
as the recent angle-averaged 3D deflagration model computed at MPA, 
Garching, by M. Reinecke et al. [29]. We believe that the main feature 
of the latter model, which allows us to get the correct flux during the 
first month, is strong mixing that moves the material enriched with 
radioactive nickel-56 to the outermost layers of the SN ejecta. 

2. Progenitors 

There is no hope to get a thermonuclear supernova from a normal star 
composed of a classical plasma. Those stars have negative effective heat 
capacity and they are thermally stable. The situation changes, if a star 
is made of a degenerate matter. The total heat capacity becomes posi- 
tive, and rTinaway can set in as in terrestrial explosives. So, a progenitor 
of SN la must be a degenerate star - a white dwarf. 

A single white dwarf is unable to explode, it cools down. But when 
it is in a binary system the chances to produce a supernova do appear 
(we need only one in ~ 300 dying white dwarfs to explode in order to 
explain the rate of SNe la). Even if the binary has two dead white 
dwarfs, it can explode because they can merge due to emission of 
gravitation waves (double-degenerate, or DD scenario [16]). If one star 
in the binary is alive, the white dwarf can accrete its lost mass and 
reach an instability (single-degenerate, or SD scenario [37, 8]). It is 
unclear which scenario is most important, there are strong arguments 
[19] from chemical evolution that only SD is the viable one. On the other 
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hand, it seems that DD can produce a richer variety of SN la events. 
Moreover, discoveries of intergalactic SNe la [4, 11] can be explained 
more naturally, because a DD system may evaporate from a galaxy. It 
is quite likely that both scenarios are being played, but their relative 
role may change in young and old galaxies. If so, a systematic trend 
may appear in SNe la properties with the age of Universe, and this 
may have important consequences for cosmology. 

3. Thermonuclear flames 

After merging in DD scenario, or after the white dwarf accretes large 

amount of material in SD case, the explosive instability develops. In 
principle, combustion can propagate either in the form of a supersonic 
detonation [1] wave, or as a subsonic deflagration [17, 24] (flame). In 
detonation, the unburned fticl is ignited by a shock front propagating 
ahead of the burning zone itself. In deflagration, the ignition is governed 
by heat and active reactant transport, i.e. by thermal conduction and 
diffusion. 

3.1. Laminar flames 

Most likely, the runaway starts as a laminar flame propagating due to 
thermal conduction. 

The rate of thermonuclear heating scales as 

(avo) ~ exp-(aG/r)^/^ 

due to the Gamow's peak: the chances to penetrate the Coulomb barrier 
for fast nuclei grow, but the tail of Maxwell distribution goes down. 
Here ac depends strongly on nuclei charges Zf ao oc ZIZ2, thus high- 
Z ions can fuse only at high T. Small perturbations of T produce huge 
variations in energy production rate since, normally, T <C aa- 
In terrestrial flames, the 'fusion' of molecules goes with the rate: 

(avo) ~ expi-Ea/TZT), 

- the Arrhenius law of chemical burning. Here Ea is activation energy. 
The parameter, showing the strong T-dependence of the heating 

Ze = 91og(CTt;o)/91ogT ~ Ea/TZT 

is called the Zeldovich number in the theory of chemical flames. For 
them typically Ze ~ 10 . . . 20. The classical theory [41] predicts the 
flame speed 

Vi^Ze-\vTlT/TreaciT2p^ 
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Table I. Flame speed Vf and width k in C+O [36] 



p 


Vi 


k 


Ap/p 


Ma 


10^ gcc 


km/s 


cm 






6 


214 


1.8 X lO"'^ 


0.10 


2 X 10-2 


1 


36 


2.9 X 10"" 


0.19 


4 X 10"^ 


0.1 


2.3 


2.7 X 10"^ 


0.43 


4 X 10"* 



with T2 - the temperature of burnt matter (ashes) and Treac(^) oc 
exp[£;„/(7^^)]. in SNe, for nuclear flames, T,eac{T) oc exp[aj^ /{3T^/^)], 
and, 

Ze = aiog((7t;o)/aiogr ~ aJ/Vl^r^^^), 

which has values very similar to terrestrial chemical flames. 

A big difference with chemical flames is the ratio of heat conduction 
and mass diffusion, the Lewis number, Le = {vtIt)/ (vdId)- Le ~ 1 
in laboratory gaseous flames, while Le ~ 10^ in thermonuclear SNe, 
since heat is transported by relativistic electrons, vt ~ c, and there is 
almost no diffusion. It ^» Id- Nevertheless, the modern computations 
[36] follow the old theory [41] closely. The conductive flame propagates 
in a presupernova with Vf which is too slow to produce an energetic 
explosion: the ratio of to sound speed, i.e. the Mach number. Ma, is 
very small (see Table I). The star has enough time to expand, to cool 
down, and the burning dies completely. So an acceleration of the flame 
is necessary in order to explain the SN phenomenon. This is the main 
problem in current research of SNe la hydrodynamics. 

3.2. Flame Instabilities and its Acceleration 

There is a rich variety of instabilities that can severely distort the shape 
of a laminar flame. The Rayleigh-Taylor (RT) instability governs the 
corrugation of the front on the largest scales. On the smallest scales the 
flame is controlled by the Landau-Darrieus (LD) instability. RT, LD 
instabilities and turbulence make computations difficult, but without 
them a star would not explode. All these instabilities were considered 
already by L. Landau [21] as a means to accelerate the flame. 

Although he did not used a term 'Rayleigh-Taylor instability'. Lan- 
dau derived his dispersion relation with account of gravity, so RT is 
there! Let k be wave number, and frequency of perturbations on the 
flame front, the shape of perturbations being of the form exp(ifey-|-r2i). 
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Here we reproduce most important expressions for the Landau disper- 
sion relation when g ^ and the surface tension of the fuel is a. The 
notation below follows Landau [21], so vi is the velocity of fuel, V2 the 
same for ashes, j = pivi = P2V2 is the mass flux across the front. Then 
the Landau dispersion relation is 

^^(^'i + ^^2) + 2^}kviV2 + k'^viV2{vi — V2) — kg{vi — V2) + ak^viV2/j = 0. 

Note that this relation contains RT instability: when k is small and Vi 
is small, we have 

Q^{vi + V2) - kg{vi - V2) = 0. 

That is, if we replace Vi by 1/pi using j = piVi = const, we get a 'pure' 
RT dispersion equation. 

When surface tension is not negligible, we have its stabilizing effect 
in RT: 

r2^(fi + V2) - kg{vi - V2) + ak^viV2/j = 0. 

But this is normally important only for very short waves, and if we 
have 

kviV2 » g , 

i.e. a case of short, but not extremely short waves, then LD becomes 
important while the term with a can still remain negligible. The latter 
inequality means that gravity can be ignored, but it is just 

i « ^ . - 

k g 

A fast flame runs from RT instability away, and the RT instability has 
no time to develop. But this does not mean that the flame becomes 
stable! No stabilization is obtained for waves below ^min, if the role of 
a or a thermal diffusion is negligible, but the instability looks now not 
as RT, but as a 'pure' LD. Since e = vi/v2 < 1 (from mass flux j = 
piVi = P2V2 and pi > P2) this gives a growing mode of LD instability: 



J^LD = --; kv2. 

e + 1 

Of course, in real life one has to use the full Landau relation, since the 
subdivision into domains of RT and LD instabilities is only approxi- 
mate. In many laboratory experiments laminar flames remain stable. 
Below a certain scale, so-called Markstein length, dissipative effects 
(like thermal conduction) dominate and smear out all perturbations, 
so that the flame may remain stable on small scales. In the majority of 
experiments the flame is attached to a burner, which sets bounds to the 
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development of LD instability at the Markstein scale, but some special 
experiments with spherical flames do show a growth of LD instability 
and fractalization of the front. Very nice examples of the action of the 
LD instability are also observed in growth of bubbles in a overheated 
liquid. 

Because of instabilities, the flame surface becomes wrinkled and its 
area grows as S cc R" , with average radius R and a > 2, i.e. faster than 
S oc R^. In other words the surface becomes 'fractal'. The exponent a 
is actually the fractal dimension, a = Dp. The effective flame speed is 
determined [39] by the ratio of the maximum scale of the instability to 
the minimum one: VeS = (AInax/Amin)■^^~^• 
Proper understanding of the behavior of flames cannot be reached if 
we restrict ourselves to a linear analysis of the problem, even if we take 
into account dissipation. However, good progress can be made with the 
help of a simple nonlinear model for the flame. 

When the vorticity is not important it is possible to study in detail 
the non-linear stage of LD instability and to find the fractal dimension 
[5]. Instead of e we define 

V = ^P/P = 1 - Ph/pu = 1 - e . 

For pi > 10^ g/cm^, which is typical for preSN cores, rj < 0.4, since 
the unburned matter is very degenerate and its thermal expansion 
relatively low. For small rj, it is easy to show that the flow behind 
the flame front remains irrotational (to first order relative to rj) if the 
incoming flow is potential and to use a simple model to compute the 
flame propagation. Laboratory flames in gases have typical values of 
rj 0.8. In this case a lot of vorticity is generated on the front and the 
simple model is not applicable. 

One can see that for low r] fractalization is not so pronounced as for 
high density jumps. We predict [5], 

Z)f(?7)(LD) ^ 1 + Z)or/^ when7?^0. 

Making a least-square fit to the power law dependence of the length of 
the fiame on the mean radius yields the values in the Table II. We find 
that Do ~ 0.3. 

3.3. Numerical Thermonuclear Flames 

The fractal description is good for LD instability while it remains mild, 
because it operates in a star on the scales from the fiame thickness (a 
tiny fraction of a cm) up to ~ 1 km. For the RT instability, Amax/ Amin is 
very uncertain and the fractal dimension is uncertain too. (Although a 
dependence of the fiame fractal dimension on the density jump across 
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Table II. Df for small 



V 




Df 


0.3 


0.09 


1.022 


0.35 


0.122 


1.039 


0.4 


0.16 


1.046 



the front was found in SPH simulations of the flame subject to RT 
instability [7] similar to the LD case). 

So a direct 3D numerical simulation is necessary. The same is true 
for a low density regime of LD when it is strongly coupled to turbu- 
lence (generated on the front, or cascading from large RT vortices). A 
great progress is achieved here in several groups [14, 29, 31, 18]. When 
simulating 3D turbulent deflagrations one encounters two problems: 
the representation of the thin moving surface separating hot and cold 
material, and the prescription of the local velocity V{ of this surface as 
a function of the large-scale flow with a crude numerical resolution > 1 
km. One solution is sketched in [29]; for a different approach sec [18]. 
In spite of the progress this problem cannot be treated as completely 
solved, and even ID approach may give interesting results, especially 
for unusual SNe la [10]. 




Figure 1. Snapshots of the flame front for 3D simulations [30]. Left panel: initial 
multi-point ignition. Right panel: flame after 0.5 seconds 
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4. Light curves of SNe la 

4.1. Basics 

At the moment, there are many models of thermonuclear explosion of 
a star, that lead to the event we know as SN la. Only a few parame- 
ters, such as kinetic energy and total ^^Ni production, can be derived 
directly from the modelling of the explosion and compared with the 
observational values. The subsequent evolution of the exploded star 
gives us much more possibilities to compare models and to decide which 
one fits observations better by reproducing more details in SN la light 
curves and spectra. We will focus here on the broad-band UBVI and 
bolometric light curve computations for SN la models. 

There are several effects in SNe physics which lead to difficulties in 
the light curve modelling of any type of SNe. For instance, an account 
should be taken correctly for deposition of gamma photons produced 
in decays of radioactive isotopes, mostly ^^Ni and ^^Co. After being 
emitted, gamma photons travel through the ejecta and can finish up in 
either thermalization or in non-coherent scattering processes. To find 
this one has to solve the transfer equation for gamma photons together 
with hydrodynamical equations. Full system of equations should in- 
volve also radiative transfer equations in the range from soft X-rays to 
infrared for the expanding medium [6]. 

There are millions of spectral lines that form SN spectra, and it 
is not a trivial problem to find a convenient way how to treat them 
even in the static case. The expansion makes the problem much more 
difficult to solve: hundreds or even thousands of lines give their input 
into emission and absorption at each frequency. 

On the first glance, the modelling of SNe la seems easier than that of 
other types of supernovae, since the hydro part is very simple: coasting 
stage starts very early, there are no shocks and no additional heating 
from them. 

On the other hand, much more difficulties arise in the radiation 
part. SNe la becomes almost transparent in continuum at the age of a 
few weeks. This means that NLTE effects are stronger than for other 
types of supernovae. Radiation is decoupled from matter within the 
entire SN la ejecta even before maximum light (which occurs around 
the 20th day after explosion), see Fig. 2 and also [12]. In this case 
one cannot ascribe the gas temperature, or any other temperature, to 
radiation, since SN la spectrum differs strongly from a blackbody one. 
Instead, one has to solve a system of time-dependent transfer equations 
in many energy groups, with an accurate prescription for treatment of 



blimisprt2.tex; 2/02/2008; 2:08; p. 8 



Type la Supernova models 



9 



0.8 - 



0.6 - 



0.2 



0.4 








0.5 



M/M, 



Figure 2. The distribution of the ratio of actual radiation energy density to the 
blackbody one at the local gas temperature inside the ejecta at different times. It 
demonstrates how radiation is decoupled from matter. 

a huge number of spectral lines, which are the main source of opacity 
in this type of SNe [3, 27]. 

Recently, powerful codes appear aimed to attack a full 3D time- 
dependent problem of SN la light [15]. Yet there are some basic ques- 
tions, like averaging the line opacity in expanding media, that remain 
controversial. 

We present theoretical UBVI- and bolometric light curves of SNe la 
for several explosion models, computed with our multi-group radiation 
hydro code. We employ our new corrected treatment for line opacity in 
the expanding medium. The results are compared with observed light 
curves. 

4.2. Method 

We compute broad-band UBVI and bolometric light curves of SNe la 
with a multi-energy radiation hydro code STELLA. Time-dependent 
equations for the angular moments of intensity in fixed frequency bins 
are coupled to Lagrangean hydro equations and solved implicitly [6]. 
Thus, we have no need to ascribe any temperature to the radiation: the 
photon energy distribution may be quite arbitrary. 

While radiation is nonequilibrium in our approximation, ionization 
and atomic level populations are assumed to be in LTE. The effect of 
line opacity is treated as an expansion opacity according to Eastman 
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t, days 



Figure 3. UBVI light curves of the model W7 with two opacity approximations: 
Eastman-Pinto (dashed) and our new approach (solid). UBV bands remain almost 
unaffected, while I-band shows now a structure more similar to what one can observe 
for real SNe la, though not identical. Crosses, stars and triangles show the light 
curves for three observed SNe la. 



& Pinto [13] and to our new recipes [33]. We have compared the results 
and found that infrared bandpass (in which the ejecta are most trans- 
parent) is more sensitive to the treatment of opacity than UBV, so one 
should be very careful on this point. See Fig. 3 for the comparison of 
SN la light curves calculated with these two approaches. To simulate 
NLTE effects we used the approximation of the absorptive opacity in 
spectral lines. NLTE results [3] and ETLA approach [27] demonstrate 
that fully absorptive lines gives us an acceptable approximation of 
NLTE effects. 

We treat gamma-ray opacity as a pure absorptive one, and solve the 
7-ray transfer equation in a one-group approximation following [35]. 
The heating by decays ^^Ni — ^ ^^Co— > ^^Fe is taken into account. 

In the calculations of SN la light curves we use up to 200 frequency 
bins and up to ^ 400 zones as a Lagrangean coordinate. 
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4.3. Light curves for presupernova models 

In our previous analysis we have studied two Chandrasekhar-mass mod- 
els: the classical deflagration model W7 [25] and the delayed detonation 
one DD4 [40], as well as two sub-Chandrasekhar-mass models: helium 
detonation model LA4 [23] and low-mass detonation model WD065 
with low ^^Ni production [32], which was constructed for modelling 
subluminous SNe la events, such as SN 1991bg. All those models were 
simplified spherically-symmetrical (ID) ones. 




T — I — I — I — I — I — I — I — I — I — I — I I — I — I — I — I — I — I — I — I — I — I — r 




t, days 

Figure 4. UBVI light curves for 4 ID models (W7 - solid line, DD4 - dots, LA4 - 
dashes, WD065 - dash-dots). Crosses, stars and triangles show the light curves for 
three observed SNe la. 

The UBVI light curves of ID models are shown in Fig. 4. The 
Chandrasekhar-mass models demonstrate almost identical light curves. 
The sub-Chandrasekhar-mass ones are more different. WD065 has al- 
most similar element distribution as Chandrasekhar-mass models, and 
the shape of its light curve is in principle the same as that of W7 and 
DD4. It is just much dimmer due to an order of magnitude less ^^Ni 
abundance. 

LA4 is very different from any other model, since the explosion there 
started on the surface of a white dwarf, not in the center, as for every 



bliimsprt2.tex; 2/02/2008; 2:08; p. 11 



12 S. Blinnikov and E. Sorokina 

Table III. Parameters of SN la models 



Model DD4 W7 LA4 WD065 MR 

Mwd'' 1.3861 1.3775 0.8678 0.6500 1.4 

M56Ni'' 0.63 0.60 0.47 0.05 0.42 

Ei,!^ 1.23 1.20 1.15 0.56 0.46 

''in Mq 

Hn 10^^ ergs s"^ 



other model, so there is a °°Ni layer near the surface in LA4. This 
feature explains why the model is essentially bluer than all three others. 
During first weeks UV quanta are locked inside the models with ^^Ni 
in the center and come out only if they are splitted into several redder 
photons, while in LA4 model they can travel outside from the surface 
freely. Unfortunately, LA4 seems much bluer than real SNe la, its rise 
and decline rates are also too fast, so centrally ignited models look 
better from the observational point of view. 

Recently, an interesting and much more sophisticated 3D deflagra- 
tion model by M.Reinecke et al. [29] (MR, hereafter) has appeared. 
Working in more dimensions the authors have to involve less free pa- 
rameters than was needed in ID simulations, so they get their model 
almost from the "first principles" . The main features of the 3D model 
are compared with the ones of classical ID models in the Table III. 

From the first glance it seems that the light curves for the models 
which differ so much could not be similar. Nevertheless, Fig. 5 demon- 
strates that they are similar in many details. The possible reasons for 
this can be understood if one has a look at the element distribution 
over the ejecta. The compositions for W7 and MR models are shown 
in Fig. 6. 

At the moment of our light curve computation the full nuclide yields 
for MR were not yet obtained. Therefore, the model consisted of the 
elements, which were chosen as representative examples for the energy 
release calculation, namely, "Fe" for iron group elements that were 
divided onto 80% of ^^Ni and 20% of ^^Fe, "Mg" for intermediate mass 
elements, and unburned C and O in equal proportion. 

The instabilities that have developed in the 3D model were not 
supposed to be so huge in approximate ID models of explosion. This 
has led to the differences in the nickel distribution over the ejecta: it is 
mixed to the outermost layers in the 3D model. These layers become 
much more opaque than in the ID models, and, despite having less 
than a half of kinetic energy, the 3D model has a photospheric velocity 
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comparable to that in the ID models. It is probably still a bit too low, 
so the light curve is wider, since the ejecta expand a bit slower, and 
photons are locked inside them for a bit longer time. The broad-band 
light curves for MR model fit the observations of one of the typical 
SN la, SN 1994D, in U and B bands surprisingly well, while classical 
ID models, such as W7 and DD4, show faster decline in the optics than 
it is observed. 

Unfortunately, the bolometric light curve for MR model (Fig. 7) is 
somewhat too slow. The ejecta must expand with higher speed to let 
photons to diffuse out faster. This means that the total energy must be 
larger to fit the bolometric observations. The bolometric light curve of 
W7 is also shown in the Fig. 7. It demonstrate much more acceptable 
decline rate, though the behaviour of the light curve before maximum 
light seems better in MR model. 

The model MR is not final. The work on getting a new, more ener- 
getic 3D model is in progress at the MPA supernova group. It seems 
that such a model should still be as much mixed as MR. Then one 
could expect that it would fit the observations of SN la bolometric 
light curves as well as the broad band ones. 
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Figure 6. Element distributions for DD4 (top left) W7 (bottom left) and MR (right) 
models. 




Figure 7. Total (thin dotted) and UVOIR (thick solid) bolometric luminosities of 
MR (left panel) and W7 (right panel) models compared with observations [9]. 



There is also another reason which allows us to believe that the MR 
model is better than the classical ID models. We calculated in details 
the X-ray emission of Tycho SNR, which is believed to be the remnant 
of SN la. The code we use takes into account the time-dependent ioniza- 
tion and recombination. We have compared the computed X-ray spectra 
and images in narrow filter bands with XMM-Newton observations of 
the Tycho. Our preliminary results [20] show that all Chandrasekhar- 
mass models produce similar X-ray spectra at the age of Tycho, but 
they differ strongly in the predictions of how the remnant should look 
like in the lines of different ions due to very different distribution of 
elements in the ejecta for ID and 3D models. We found that W7 and 
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DD4 models produce rather wide ring in Fe lines, while it is narrow for 
MR model. The image for the latter model is very similar to what is 
observed. 

We believe that the main feature of this model which allows us to 
get correct radiation during the first month, as well as after a few hun- 
dred years, when an SNR forms, is strong mixing that pushes material 
enriched in iron and nickel to the outermost layers of SN ejecta. 

5. Conclusions 

There are many points which require attention in research of SNe la: 
progenitors; burning regimes (that may change with the age of Universe 
[34]). The physical understanding of the Pskovskii-Phillips is not yet 
achieved (probably it will be reached when the burning will be modeled 
completely from the first principles, because too many parameters enter 
in light curve computations). 

The new 3D SN la model MR [29] is very appealing. Yet it is not 
a final one: a detailed post-processing of nucleosynthesis changes the 
composition. It has been done very recently [38], and it is not yet 
checked in the light curve calculation. Our light curve computations are 
also preliminary, since more work is needed on the expansion opacity. 
Hopefully, none of the required improvements will spoil the light curve 
of this model and its X-ray spectra on the SNR stage, since the specific 
qualities of the model can be primarily explained by the enrichment of 
the outermost layers of SN ejecta by Fe and Ni. 

The SN light curve modelling still has a lot of physics to be added, 
such as a 3D timc-dcpcndcnt radiative transfer, inchiding as much as 
possible of NLTE effects, which are especially essential for SNe la [15]. 
All this will improve our understanding of thermonuclear supernovae. 
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